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ABSTRACT 

We investigate in massive early-type galaxies the variation of their two-dimensional central 
fraction of dark over total mass and dark matter density as a function of stellar mass, central 
stellar velocity dispersion, effective radius, and central surface stellar mass density. We use a sam- 
ple of approximately 1.7 x 10 5 galaxies from the Sloan Digital Sky Survey Data Release Seven 
(SDSS DR7) at redshift smaller than 0.33. We apply conservative photometric and spectroscopic 
cuts on the SDSS DR7 and the MPA/JHU value-added galaxy catalogs, to select galaxies with 
physical properties similar to those of the lenses studied in the Sloan Lens ACS (SLACS) Survey. 
The values of the galaxy stellar and total mass projected inside a cylinder of radius equal to the 
effective radius are obtained, respectively, by fitting the SDSS multicolor photometry with stellar 
population synthesis models, under the assumption of a Chabrier stellar initial mass function 
(IMF), and adopting a one-component isothermal total mass model with effective velocity dis- 
persion approximated by the central stellar velocity dispersion. The plausibility of an isothermal 
model to represent the galaxy total mass distribution is supported by independent gravitational 
lensing and stellar dynamical analyses performed in the lens subsample, which is found here to 
represent nicely the entire galaxy sample. We hnd that within the effective radius the stellar mass 
estimates differ from the total ones by only a relatively constant proportionality factor. In detail, 
we observe that the values of the projected fraction of dark over total mass and the logarithmic 
values of the central surface dark matter density (measured in Mq kpc -2 ) have almost Gaussian 
probability distribution functions, with median values of 0.64+q'5i and 9.1+q' 2 ,) respectively. We 
discuss the observed correlations between these quantities and other galaxy global parameters 
and show that our results disfavor an interpretation of the tilt of the Fundamental Plane (FP) in 
terms of differences in the galaxy dark matter content and give useful information on the possible 
variations of the galaxy stellar IMF and dark matter density profile. Finally, we provide some 
observational evidence on the likely significant contribution of dry minor mergers, feedback from 
active galactic nuclei, and/or coalescence of binary black holes on the formation and evolution of 
massive early-type galaxies. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: structure — dark matter 
1. Introduction 
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chical model of structure formation, i.e., the con- 
cordance ACDM (co l d dark matter) mode l (e.g., 
Springel et ail 120061; iKomatsu et all l2009h . De- 



spite this success, at smaller scales (« 1 kpc) a 
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complete understanding of the dark and baryonic 
matter interplay, that ultimately determines the 
internal structure of galaxies, is still missing. In- 
terestingly, not only the amount and distribution, 
but also the presence of dark matter in early-type 
galaxies has been controversial for a long time. 
The lack of suitable and easily interpreted kine- 
matical tracers, such as Hi in spirals, at large 
radii and the degeneracy between the mass dis- 
tribu tion and the anisotropy of the stellar orbits 
(e.g., lGerhardl[l995 lLokas fe Mamonll2003l) have 
made the detection of dark matter through dy- 
namical studies difficult. Spectroscopic observa- 
tions of the light emitted by stars have though 
provided the first useful constraints on the mass 
distribution within approximately two times the 
value of the eff ective radius R e in lo c al early-type 
galax i es (e.g.. | van der Maxell Il99ll: iBertin et al 



19941: ICarollo et all 1 199.4 iGerhard et alJ 1200 it 
Cappellari et all " 20061) . In the framework of 
collisionless collapse, fully self-consistent two- 
component (luminous and dark ma tter) dynam- 



ical models have been developed (|Bertin et al 



19921 ) and applied (jSaglia et al1ll992l ) to fit the 



photometric and spectroscopic data of a sam- 
ple of nearby galaxies, indicating positive evi- 
dence for the presence of dark matter. More 
recently, Schwarzschild (orbit-based) dynamical 
models have also enabled a decomposition of the 
total mass into luminous and dark in a sample of 
Coma early- type galaxies (jThomas et al.ll2007f) . 

Other kinematical tracers, like cold atomic hy- 
drogen and warm ionized gas, have supported the 
picture of a dark matter component additional 
to the luminous one in early-type ga laxies (e.g., 
iBuson et all Il993t iFranx et all Il994l ). Unfortu- 
nately, these gaseous sources are either rare or too 
restricted radially to trace dark matter straight- 
forwardly in statistically significant samples of 
galaxies. A different diagnostics of the gravita- 
tional potential is the diffuse X-ray emission from 
the hot gas present in many early-type galaxies. 
The properties of the gas have allowed in some 
cases to estimate, under the hypotheses of quasi- 
hydrostatic equilibrium and spherical symmetry, 
mass-to-light ratios on the orde r of 100 on radial 
scale s of about 100 kpc (e.g., iMushotzkv et al 



1994; lLoewenstein fc W hite 1999). Two difficul 



gas is complex to model and the temperature pro- 
file of the gas may be measured with considerable 
uncertainty. A different class of tracers includes 
globular clusters and planetary nebulae. The dis- 
crete kinematical data obtained by observing these 
objects orbiting around some nearby early-type 
galaxies have provided an estimate of the grav- 
itational field in these galaxies, co nfirming the 



presence of dark matter ha los (e.g., I Mould et al 
19901 : lArnaboldi et al.lll998l ). All the cited anal- 
yses agree reasonably well in revealing flat cir- 
cular velocity curves (i.e., a total matter density 
distribution well approximated by a 1/r 2 profile) 
from > 0.2 to > 2 R e in early-type galaxies. Not 
only does the associated mass-density profile differ 
signi ficantly from cosmol ogically motivated ones 
(e.g., iNavarro et al 1996 ). but it also requires a 
significant amount of fine-tuning (known as bulge- 
halo "conspiracy" ) between the distribution of lu- 
minous and dark matter. This high degree of ho- 
mology is still poorly understood in the currently 
accepted scenario of galaxy formation and evolu- 
tion. 

During the past decade, gravitational lensing 
has been extensively used to study the mass dis- 
tribution of early-type galaxies beyond the local 
Universe, con tributing to key observations of dark 
matter (e.g., Rusin et al. 2003: Koopmans et al 



ties affect studies of this kind, a theoretical and an 
observational one: the physics of the hot cooling 



20061 : lOriHo et alJl2009t IVegetti et al.fl2010h . The 

combination of strong lensing and stellar kine- 
matics has proved to be particularly effective, 
since the two diagnostics complement each other 
Ce.g..lTreu fe Koopmansll2004t ICzoske et alll2008t 
Barnabe et al.ll2009l ). In fact, strong lensing pro- 
vides a robust measurement of t he total mass pro - 
jected inside the Einstein radius (|Kochaneklll99ll) , 
breaking the mass-anisotropy degeneracy of the 
dynamical analysis, after which stellar dynamics, 
giving constraints on the total mass distribution at 
small radii (typically < i? e ), provides an estimate 
of the mass gradient. Thus, a joint strong lensing 
and stellar-dynamical study can be used to deter- 
mine the average logarithmic density slope of the 
total mass inside the Einstein radius of individual 
lens galaxies, independent of the mass- sheet de - 
generacy (e.g., Falco et al. 1985 : Kochanek 20061 ). 
In addition, this allows one to decompose the 
total mass distribution into luminous and dark 
with good precision, making feasible an investi- 
gation on the evolution of the internal structure 
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of early-type galaxies. This technique has been 
applied to a sample of 5 massive early-type galax- 
ies in the redshift range z « 0.5 — 1, as part of 
the Lenses Structure and Dynamics (LSD) Survey 
( Treu fc Koopmansl 2004 ). and to a larger sam- 
ple of 58 massive early-type galaxies at lower red- 
shifts (z w 0.06 - 0.5) from the Sloan Le ns ACS 



invest igation of the most plaus i ble stellar IMF 



(SLACS) Survey (|Koopmans et al.1 120091) . De 



spite some scatter in individual galaxies, one of 
the most important results of these two surveys is 
the measurement in the complete sample of a re- 
markably homogeneous total (luminous and dark) 
mass density profile, that is consistent with an 
isothermal one (i.e., p oc 1/r 2 ) out to a few hun- 
dreds of kiloparsecs (see also lGavazzi et al.ll2007 



Bolton et al. 12008b ). Additional support to the 



picture of an average isothermal total mass den- 
sity profile with a surprisingly small scatter has 
been provided by a fully self-consistent analysis of 
lensing data and two-dimensio nal kinematic maps 
in a smaller subsam ple (e.g., Czoske et al.l [2008; 
Barn abe et al. l2009h . Moreover, no evidence of 



significant evolution in the value of the total mass 
density slope below redshift of 1 has been found. 
The measured one-component isothermal model is 
completely characterized by the value of an ef- 
fective velocity dispersion, which has resulted to 
be nicely approximated (within less than 3%) by 
the value of the galaxy central stellar velocity dis- 
persion (i.e., the projected velocity dispersion of 
the stars within a disk of radius R e / 8; see, e.g., 
Treu et al.l 120061: iBolton et ail l2008bh . This lat- 
ter result is expected when considering the Jeans 
equat ion for real is tic stellar densit y distributions 
(e.g., IJaffd 119831 iHernquisd fl99oT ) em bedded in 
a glo bally isothermal distribution (see iKochanek 
1993) and some observational evidence has been 
collected by means of dynamical modeling in sam- 
ples of nearby and luminous early-type galaxies 
fe.g.. lKochanek1ll994tlGrillo et al.ll2008bh . 

Strong lensing has also turned out to be an in- 
valuable astrophysical tool when combined with 
stellar population synthesis models. Taking ad- 
vantage of this combination, recent studies have 
addressed disparate topics that are relevant in the 
field of galaxy formation and evolution, such as the 
measurement of the projected dark over total mass 
fractions inside the i nner regions of dista nt early- 



type galaxies Ce.g IGrillo et al.l l2008allcl. J2009l 



2010al : lAuger et all 120091: iFadelv et al.1120101 ). the 



(e.g., IGrillo et all l2008al 120091 : IGrillo fc Gobat 



2010b 



rcu ■ 



_ et al.l 2OIOI) and origin of the tilt o f 
the FP (IGrillo et al.ll2009l : IGrillo fc Gobatll2010bh 

of massive early- type galaxies. Moreover, dif- 
ferent analyses have revealed that the SLACS 
lens galaxies are a representative sample of their 
SDSS parent sample with regard to photomet- 
ric, spectroscopic, and enyironmental properties 



(Bolton et al. 2006, 2008a; Treu et al.ll200 



operties 
6l l2009l 



Grillo et al.l 20091) 7 This last point encourages to 



adopt here a different perspective and exploit, in 
the light of the latest lensing results, the wealth of 
information already accumulated by the SDSS to 
study the projected central dark over total mass 
fractions and dark matter densities of massive 
(lens and non lens) early-type galaxies. 

The paper is organized as follows. In Sect. 
2, we describe how the sample of massive early- 
type galaxies is selected from the SDSS DR7 and 
the relevant galaxy photometric and spectroscopic 
values are estimated. In Sect. 3, we outline the 
method used to measure the central dark over to- 
tal mass fractions and dark matter densities. In 
Sect. 4, we present our results and explore pos- 
sible correlations between the dark matter and 
several other physical quantities of the galaxies 
in the sample. In Sect. 5, we discuss and in- 
terpret the results of our analysis within differ- 
ent galaxy formation and evolution scenarios. Fi- 
nally, in Sect. 6, we summarize and draw con- 
clusions. The aperture mass values discussed in 
this work are intended two-dimensional mass mea- 
surements, i.e., projected along the line of sight. 
All length, mass, velocity, and surface mass den- 
sity values considered in our study are expressed, 
respectively, in units of kpc, M Q , km s _1 , and 
M©kpc~ 2 . The logarithms assume a base of 10 
and dimensionless arguments obtained by divid- 
ing the studied quantities with their correspond- 
ing meas urement unit s. Following the general so- 
lution of [York! (|1966l ) to the linear least-squares 
problem, we determine all the best-fit correlation 
lines between two variables taking into account 
the uncertainties in both coordinates. Through- 
out this work, we assume H — 70 km s^ 1 Mpc -1 , 
fl m = 0.3, and fl A = 0.7. 
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Fig. 1. — The probability distribution functions of 
the logarithms of the total stellar mass M* (on the 
left) and central stellar velocity dispersion cr (on 
the right) of the approximately 1.7x 10 5 galaxies in 
the sample. The stellar masses and stellar veloc- 
ity dispersions have been estimated, respectively, 
from the SDSS photometric and spectroscopic ob- 
servations. 



2. The sample 

The sample of galaxies analyzed in this paper 
is drawn from the SDSS. This survey is providing 
multi-b and photometry in fiv e bandpasses (u, g 1 r, 
i. z: see 



Fukugita et al.lll996l ) for almost a quarter 



of the sky at high Galactic latitude and spectro- 
scopic obs ervations (taken usin g 3-arcsec diameter 



fibers; see lStrauss et al.l 120021 ) for more than 10 
objects. In addition to the data retrieved from the 
SDSS DR7 Galaxy catalog, we use here also the 
public galaxy catalogs provided by the MPA/JHU 
collaboration^. 

First, we match the catalogs and select early- 
type galaxies by requiring that the surface bright- 
ness profile of each object is well described by a de 
Vaucouleurs profile. In detail, we make use of the 
SDSS morphological index fracDeV, which quan- 
tifies the weight of a de Vaucouleurs profile in a 
two-component (a de Vaucouleurs plus an expo- 
nential profile) decomposition of the galaxy sur- 
face brightness distribution. We apply a highly 
conservative cut choosing only those objects that 
exhibit r, i, and z fracDeV values greater than 
0.95. Moreover, to make the sample more homoge- 
neous we exclude the galaxies that belong to the 
subclasses called AGN, BROADLINE, and STAR- 
FORMING in the MPA/JHU galaxy catalog. 

Then, we decide to concentrate our study 



only on massive early-type galaxies with phys- 
ical scales similar to those of the lens galaxies 
studied in the SLACS survey (see Bolton et al.l 
2008at iGrillo et al.l 120091: 1 Auger et aJll2009h . Ac- 
cordingly, we select the galaxies that have SDSS 
spectroscopic redshifts z sp between 0.05 and 0.5, 
aperture stellar velocity dispersions <7 ap that range 
from 150 and 400 km s _1 , and total stellar masses 
M* between 10 10 5 and 10 12 M Q . Stellar masses of 
the SDSS galaxies are measured by the MPA/JHU 
team by fitting the SDSS broad-band photom- 
etry (modelMag magnitudes) with a large grid 
of stellar population synthesis models, that can 
accommodate various st ar formation histories. 
Bruzual & Chariot (2003) templates are used and 
a IChabrierl (|2003l ) stellar IMF is assumed (for 
mo re details on the stellar mass measurem ents , 
see iKauffmann et alJEool ISalim et alJl2007t ). In 
the catalog, the stellar mass estimates are avail- 
able only for galaxies at redshift smaller than 0.33. 
Therefore, we restrict our sample to this upper 
limit in redshift. 

We compute the rest frame r-band effective an- 
gles 9 e (i.e., the angles within which half of the 
total light of a de Vaucouleurs profile is included) 
of the galaxies in the sample by linear interpo- 
lation in wavelength between the effective angles 
e ,r, &e.i, and 6 e . z that are fitted, considering the 
appropriate values of the point spread function, on 
the galaxy surface brightness distributions in the 
r-, i-, and z-band, respectively. For these bands, 
we take their central wavelengths A r , A^, and X z 
to be correspondingly 6231, 7625, and 9134 A, re- 
sulting in these expressions for the galaxy angular 
sizes: 



= e e . r + [(i + z sp )\ r - x r ] ■ 



Xi-Xr 



1 http: / /www. mpa-garching.mpg.de/SDSS/ 



0e = 6e,i + [(1 + 2sp)A r - A,-] %f^f- (1) 

The first and second formulae shown in Eq. ([T]) 
are used, respectively, if the spectroscopic redshift 
of a galaxy is smaller or larger than 0.224. The 
effective radii R e are determined by multiplying 
the values of the effective angles by the angular 
diameter distances at the galaxy redshifts. 

Next, we estimate the central velocity disper- 
sions of the galaxy stellar component, coj start- 
ing from the SDSS velocity dispersions, cr ap , that 
are measured within a 3-arcsec diameter fiber, and 
using the following empirical prescription deter- 
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mined by I J0rgensen et aL ( 19951) 



8 x 1.5" 



-0.04 



(2) 



The sample of massive early-type galaxies that 
we select according to the selection criteria de- 
scribed above consists of approximately 1.7 x 10 5 
objects. In Fig. [TJ we plot the resulting probabil- 
ity distribution functions of the logarithms of the 
mean values of the total stellar mass and of the 
values of the central stellar velocity dispersion. 

We emphasize that the total stellar mass esti- 
mates used in this study are robust. In fact, these 
quantities, that are measured from the SDSS mul- 
ticolor photometry, agree very well (especially at 
the high mass end considered here) with the values 
determined by fitting different spectral indices in 
the galaxy spectre ^. Moreover, in a recent study 
Auger et all |2009) have compared in a subsample 
of approximately 50 massive lens galaxies the pho- 
tometric stellar masses estimated by modeling in 
different ways (as done by the M PA/JHU collabo- 
ration and bv lGrillo etal1l2009l) the public SDSS 
photometry and their new HST photometry. De- 
spite some minor differences, the authors conclude 
on the overall consistency of all the independently 
measured stellar masses. 

3. The method 

We focus here on the determination of two- 
dimensional dark over total mass fractions and 
dark matter densities, measured within cylinders 
of radii R e . We use the luminous (stellar) mass 
values introduced in the previous section and to- 
tal mass values estimated from stellar dynamics. 
We intentionally concentrate on projected quanti- 
ties for two reasons: first, because these quantities 
are more directly related to the observables, with 
no need of any deprojection modeling, and second, 
because these are the physical quantities that en- 
ter in several scaling laws of early-type galaxies, 
like the FP. 

Based on the assumption that the light distri- 
bution of the galaxies traces well their stellar mass 



2 A discussion on the observed differences between 
the estimates inferred from the photometric and 
spectroscopic stellar mass diagnostics is given at 
http: / /www.mpa-garching. mpg.de/SDSS/DR7/mass_comp. html 



distribution and on the approximation that half of 
the galaxy total light is enclosed within a disk of 
radius equal to R e , we estimate the values of the 
stellar mass within the effective radius M*(< R e ) 
by halving the photometrically determined mean 
values of the total stellar mass M*: 



M*(< R e ) = M./2. 



(3) 



The stellar masses are known with typical errors 
on the order of 20%. 

Then, assuming that the total (luminous and 
dark) mass distribution is well described by a one- 
component isothermal distribution and that the 
value of (To is representative of the value of the ef- 
fective velocity dispersion of the isothermal model 
(for more details on the observational evidence rel- 
ative to these points, see Sect. 1 and references 
therein) , we measure the values of the total mass 
within the effective radius Mt(< R e ) as follows: 



M T (< Re) = Tra 2 R e /G. 



(4) 



The total masses have a typical error on the order 
of 12%. 

Next, we compute inside the effective radius the 
fraction of dark over total mass /d(< -Re) by sub- 
tracting to one the luminous over total mass frac- 
tion inside R e : 

/ D (< Re) = 1 - M»(< R e ) I M T (< Re) , (5) 

and the central surface dark matter density £d,-R c 
by dividing the estimate of the mass in the form 
of dark matter by the area of the disk delimited 
by R e : 

£da = [M T {< Re) - M*{< Re)] /nR 2 e . (6) 
4. Results 

We plot in Fig. Hthe 68%, 95%, and 99% confi- 
dence regions of the luminous versus total masses 
and the dark matter fractions versus dark matter 
densities (all the quantities are estimated within 
i? e ), and the probability distribution functions of 
the dark over total mass fractions and dark mat- 
ter densities of all the galaxies in the sample. In 
each panel of Fig. [2] we also show the 44 lens 
galaxies of the SLACS survey that according to 
the adopted selection criteria belong to the sam- 
ple studied here. 
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Fig. 2. — On the top left: The logarithmic values of the galaxy total Mr(< R e ) and stellar M*(< R e ) 
masses projected in a cylinder with radius equal to the effective radius R e . On the top right: The probability 
distribution function of the fractions of dark over total mass projected within the effective radius /d(< Re)- 
On the bottom left: The probability distribution function of the logarithmic values of the dark matter density 
projected within the effective radius Sd,k c • On the bottom right: The values of the fraction of dark over 
total mass and the logarithmic values of the dark matter density projected in a cylinder with radius equal to 
the effective radius R e . For the approximately 1.7 x 10 5 galaxies in the sample, the contour levels represent 
the 68%, 95%, and 99% confidence level regions, and the solid and dashed lines show the best-fit and the 
one-to-one lines, respectively. The solid lines on the i-axis represents the 68% confidence intervals. The 
squares and the light color probability distribution functions indicate the 44 (grade-A) SLACS lens galaxies 
that satisfy the selection criteria of the sample. 
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As expected, we find a highly statistically sig- 
nificant value of 0.91 for the Pearson linear cor- 
relation coefficient between the logarithmic values 
of M T (< R e ) and M*(< R e ). The best-fit corre- 
lation line between the logarithmic values of the 
luminous and total masses of the galaxies in the 
sample is 

log [M T (< R e )] = -0.58 + 1.09 X log [M*(< Re)]. 

(7) 

The errors on the two best-fit linear coefficients are 
very small, given the large number of data points 
and their relatively small scatter around the best- 
fit line. Remarkably, if we compare the best-fit and 
one-to-one lines in Fig. [5J we can conclude that 
inside R e the galaxies in the sample show total 
masses that increase almost linearly with the lu- 
minous masses. This result is emphasized by look- 
ing at the second panel of Fig. [5] The dark over 
total mass fractions present a unimodal, slightly 
asymmetric, and almost Gaussian probability dis- 
tribution, with a median value of 0.64 and 68% 
confidence level values between 0.53 and 0.72. The 
logarithmic values of the dark matter density also 
exhibit a probability distribution function which 
is well approximated by a Gaussian function, with 
a median value of 9.1 and a 68% confidence level 
interval extending from 8.9 to 9.3. We measure a 
Pearson linear correlation coefficient value of 0.49 
between the values of the dark matter fraction and 
the logarithmic values of the dark matter density 
(for the best-fit line parameters, see Table [T] and 
Fig. [2J. This correlation suggests that an increas- 
ing fraction of dark over total mass within R e is as- 
sociated to an increasing concentration inside the 
same radius of the dark matter profile. 

The tight selection criteria adopted to define 
our galaxy sample guarantee the robustness of 
the previous conclusions. In fact, the requirement 
on the good approximation of the galaxy surface 
brightness distribution by a de Vaucouleurs pro- 
file ensures that the effective radius is suitable to 
represent the half-light radius. As a result, the 
value of the stellar mass within R e is well repre- 
sented by half the value of the total stellar mass. 
In addition, the relatively small ranges of allowed 
values for the stellar velocity dispersion and stellar 
mass make certain that a one-component isother- 
mal distribution is an appropriate parametrization 
of the total mass distribution. This last statement 
is supported by the results of numerous dynamical 



and lensing studies in smaller samples of galaxies 
that present structural properties similar to those 
of the galaxies in our sample (see Sect. 1). 

To further check how the previous result of an 
almost constant dark over total mass fraction is 
affected by some scatter in the one-to-one relation 
between the effective velocity dispersions of the 
one-component isothermal model and the central 
stellar velocity dispersions, we perform the follow- 
ing test. We extract new values of Co from Gaus- 
sian probability distributions centered on the mea- 
sured values and with standard deviations equal to 
10% of the measured values. Then, we use them to 
estimate the projected total masses and dark over 
total mass fractions, as described in the previous 
section. We find that the approximately linear re- 
lation observed in the first panel of Fig. [2]between 
the logarithmic values of the luminous and total 
masses of the galaxies in the sample is conserved, 
despite some expected larger scatter. The proba- 
bility distribution function of the dark over total 
mass fractions plotted in the second panel of Fig. 
[2] is slightly broader, but not significantly modi- 
fied. The median value is still 0.64 and the 68% 
confidence level interval extends between 0.49 and 
0.74. On the basis of this analysis, we confirm 
the robustness of our results. We remark that 
this test also ensures that small deviations from 
an average one-component isothermal model for 
the galaxy total mass distribution would act in the 
same way and not change appreciably the general 
picture presented so far. 

We emphasize that the uncertainties on the to- 
tal and stellar masses are not correlated. This en- 
sures that the correlation pattern visible in Fig. [2] 
between Mt(< R e ) and M*(< R e ) is not affected 
by the error correlations. In fact, on the one hand, 
we have checked that the errors on the stellar 
masses are primarily attributable to the degen- 
eracies that are inherent in the stellar population 
modeling and only to a secondary level to the pho- 
tometric uncertainties. On the other hand, we ob- 
serve that the errors on the total masses reflect 
essentially the precision with which the stellar ve- 
locity dispersions can be measured from the galaxy 
spectra. 

In passing, we me ntion that the baryon ic Tully- 
Fisher relation (e.g. jMcGaugh et al 1 l2000h . a scal- 
ing law between the total baryonic mass and the 
maximum rotation velocity of late-type galaxies, 
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seems also to suggest that in these objects the 
dark over total mass fractions are constant over 
a wide mass range. 

Recent analyses that combine photometric, dy- 
namical, and lensing data in massive early-type 
galaxies seem to agree on favoring for these objects 
a Salp eter to a Chabrier o r a Kr o upa stellar IMF 
(e.g.. IGrillo et all l2008aL 120091: iGrillo fc Gobatl 



2010l ). According to these re 



2010bl; iTreu et al 
suits, if we normalize in our sample the stellar 
mass estimates obtained by adopting a Chabrier 
IMF to a Salpeter IMF, we measure a median 
value of the dark over total matter fraction of 
approximately 0.4 within R e . This result is con- 
sistent, given the errors, with the average val- 
ues of 0.3 (see Koopmans et al.|[2006l : IGrillo et~al 



2009: 



Auger et alj|2009h and 0.4 fsee lrBolton et al 



2008bJ) for the 2D dark over total mass fractions 
determined, respectively, inside disks with aper- 
tures equal to the average Einstein (0.6 times the 
effective radius) and effective radii of the SLACS 
lens galaxies. Furthe rmore, in 6 of these galaxies 
Barnabe et all (|2009l ) have combined integral field 
spectroscopy with lensing data and measured val- 
ues between 15% and 30% for the 3D dark over 
total mass fractions within R e . Our results are 
consistent with these values and with the values 
between 20% and 30% for the same fractions that 
have been estimated by precise dynamical mod- 
eling of the photometric and spectroscopic obser- 
vations of about 2 massive early-type galaxies of 



the Coma cluster ( Thomas et al.ll2007l ). Unfortu- 



nately, a direct comparison with th e values of the 



dark m atter densities provided by I Thomas et al 
(2009) and determined in the Coma galaxies with 
the same dynamical technique is not feasible be- 
cause of projection effects and differently probed 
spatial regions. In fact, the dark matter densities 
of the Coma sample are measured within a sphere 
of radius equal to two times the value of R e . 

In Figs. [3] and S] and Tabled! we show the rela- 
tions, the best-fit lines, the best-fit linear parame- 
ters (a and b), and the values of the Pearson linear 
correlation coefficient (p) between the dark over 
total mass fractions and the logarithmic values 
of the dark matter densities versus the logarith- 
mic values of the total stellar mass, central stel- 
lar velocity dispersion, effective radius, and sur- 
face stellar mass density within the effective radius 
[£*,fl e = M*(< R e ) /ttRI}. We observe that the 



values of /d(< -Re) are correlated at a higher sig- 
nificance level to the values of <jq and £*.,r c than 
to those of M* and R e . The values of SD,i?„ are 
correlated at a statistically significant level to the 
values of 00, Re, and and at a low level 

to those of M*. The almost linear relation be- 
tween Xd..r c and is essentially another way 
of showing that within the effective radius both lu- 
minous and dark masses differ from the total mass 
by nearly constant proportionality factors (see also 
the first panel of Fig. [I). 

We notice that in Figs. [2]|4] the subsample 
of 44 massive lens galaxies of the SLACS sur- 
vey represents well the general properties of the 
entire sample. We remark that the SLACS lens 
galaxies are to a first approximation selected ac- 
cordi ng to the values o f their velocity dispersion 



Bolton et al. 20061 ). This justifies the pres- 



(see 

ence in our plots of slightly more lenses in the 
regions where the velocity dispersions are high. 
Moreover, the probability that a galaxy produces 
multiple images of background sources increases if 
the galaxy central tota l mass density is high (e 



van de Ven et al. 2009; Mandclb aum et all 1200 



This explains the observed moderate preference of 
the lenses to lie in regions characterized by high 
central mass densities. 

Surprisingly, our previous result on a fairly con- 
stant dark over total mass fraction within the ef- 
fective radius is at variance with some other stud- 



les (e .g., iPadmanabhan et al.l 12004c iTortora et al 
2009). A first plausible explanation of this fact 



is that in these works the dark matter fractions 
of early-type galaxies are given within a three- 



Table 1: Best-fit parameters of the dark matter 
scaling relations (y = a + b x x). 



y 


X 


a 


b 


P 


/ D (< Re) 


log(M*) 


1.27 


-0.05 


-0.12 


/ D (< Re) 


log (00) 


-0.44 


0.47 


0.30 


/ D (< Re) 


l0g(#e) 


0.60 


0.11 


0.18 


/ D (< Re) 


log(£*,flJ 


2.45 


-0.20 


-0.45 


/ D (< Re) 


log(£ D ,i?J 


-1.04 


0.19 


0.46 


log(£ D ,flJ 


Iog(M«) 


11.8 


-0.24 


-0.29 


log(£ D ,R e ) 


log(cr ) 


4.38 


2.03 


0.48 


log(£ D ,R e ) 


logCRe) 


9.69 


-0.68 


-0.54 


log(£ D ,flJ 


log(£*,H e ) 


1.17 


0.90 


0.57 



8 




log(R e ) logGVO 



Fig. 3. — The values of the fraction of dark over total mass projected within the effective radius /d(< Re) 
as a function of the logarithmic values of the total stellar mass (on the top left), central stellar velocity 
dispersion <tq ( on the top right), effective radius R e (on the bottom left), and central surface stellar mass 
density £*,_R e (on the bottom right) for the approximately 1.7 x 10 5 galaxies in the sample. The contour 
levels represent the 68%, 95%, and 99% confidence level regions, and the solid lines show the best-fit lines. 
The squares indicate the 44 (grade- A) SLACS lens galaxies that satisfy the selection criteria of the sample. 
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dimensional radius and estimated by assuming 
some virial relation or using non-trivial dynam- 
ical models. Analyses of this kind are in gen- 
eral considerably more mo del- dependent than our 
study. Furthermore, the broader and more hetero- 
geneous physical properties of some galaxy sam- 
ples (often extending to significantly less massive 
galaxies) may be a possible source of confusion 
in the interpretation of the results. We have 
considered here only massive early-type galaxies 
that have been shown to be a fairly uniform class 
of galaxies in many respects, not least as far as 
their age and metall i city c ontent is concerned (e.g., 
Gallazzi et aT1l2005L liooi ). 



5. Discussion 

We reconsider here some of the relations pre- 
sented in the previous section and explore the pos- 
sible physical mechanisms that are at their origin 
and the implications in the framework of galaxy 
formation and evolution studies. 

Early-type galaxies are found to occupy only 
a small fraction of the parameter space de- 
fined by the three variables effective radius, cen- 
tral stellar velocity di spersion, and effective sur- 
face brightness (e.g ., Diorgovski fc Daviesl 1987 ; 
Dressier et al. I Il987l) . The observed tight scal- 
ing law between these global properties, known 
as the FP, can be interpreted in terms of a 
systematic increase of galaxy effective (dynam- 
ical) mass-to-lig ht ratio with effective (dynami- 
cal) mass (e.g., Faber et all 1987t Bender et al 



19921 : Ivan Albada et al.lll995HCiotti et al.lll996l ) 



This interpretation is usually referred to as the 
tilt of the FP relative to the expectation of the 
virial theorem. The origin of the tilt of the 
FP is still debated and it can be ascribed to 
variations in the stellar population properties, 



ties of early-type galaxies (e.g., Hiorth & IV 


adsen 


1995; Prugniel & Simienlll997l; Bertin et al. 


2002 




Truiillo et al. 


2004; Padmanabhan et al. 


2004 


Grillo et al.l2009llAllanson et al.ll2009l;lTortora et al. 


20091: Grillo & Gobat 2010bh. 



Our measurement of an almost constant ratio 
between the luminous and total mass inside the 
effective radius excludes that the tilt of the FP 
of massive early-type galaxies is primarily due to 
variations in their dark matter fractions. This 



fact combined with the observed small deviations 
from homology in the mass-dynami cal structure 
of massive early-type galaxie s (e.g., iBolton et al 



2008bl ICappellari et al.ll2006l ) support the conclu- 



sions of different analyses on a main stellar pop- 
ulation origin for the tilt of the FP of these ob- 



jects fe.g..lGrillo et all liool lAllans on et al.ll2009t 
Grillo fc Gobatll2010bl ). 



A few studies have claimed to have provided 
some observational evidence that the stellar IMF 
of early- type gala xies may not be universal and 
non-e volving (e.g.. Ivan Dokkuml liooi iTreu et al 
2010h . If we followed these suggestions and ex- 



plicitly allowed a change in the stellar IMF go- 
ing from a Chabrier/Kroupa-like to a Salpeter-like 
IMF, moving from less to more massive galaxies in 
our sample, we would conclude that in the inner 
regions less massive galaxies are more dark-matter 
dominated than more massive ones. This is the 
opposite of what would be required to explain the 
tilt of the FP with varying the galaxy dark matter 
content. Alternatively, if we considered a variation 
in the stellar IMF for the galaxies in our sample, 
but fixed the previous result on an almost constant 
value of /d(< Re), this would leave small place to 
only a very fine-tuned combination of variations 
in the stellar IMF and dark matter density pro- 
file. In particular, a change in the direction from 
a Chabrier/Kroupa to a Salpeter stellar IMF as 
described above, i.e., moving towards increasing 
galaxy masses, could only be well coupled with an 
increase in the concentration of the dark matter 
density distribution to preserve the dark over to- 
tal mass fraction in the central regions. 

N-body cosmological simulations have pre- 
dicted that in an expanding universe cold dark 
matter particles collapse into self-similar ha- 
los with a diverging inner density p rofile (e.g., 
Navarro et all 19961 : [Moore et al.lll999h . The more 



recent inclusion of the baryonic physics in simu- 
lations has shown that the luminous component 
can modify signifi cantly the dark matter pr ofile 
of a galaxy (e.g., lLackner fc Ostriker 2010l and 
references therein). The assembly of the stellar 
mass of an elliptical galaxy is not easy to simu- 
late and involve different physical processes that 
can be simplistically classified as dissipationless 
and dissipational, depending on if the total en- 
ergy of the system is conserved or not. Two pos- 
sible formation scenarios are one in which stars 
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Fig. 4. — The logarithmic values of the surface dark matter density projected within the effective radius 
^D,K e as a function of the logarithmic values of the total stellar mass M* (on the top left), central stellar 
velocity dispersion cto (on the top right), effective radius R e (on the bottom left), and central surface stellar 
mass density S*,_r c (on the bottom right) for the approximately 1.7 x 10 5 galaxies in the sample. The contour 
levels represent the 68%, 95%, and 99% confidence level regions, and the solid lines show the best-fit lines. 
The squares indicate the 44 (grade-A) SLACS lens galaxies that satisfy the selection criteria of the sample. 
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originally form far from the effective radius of a 
galaxy and are subsequently accreted in its inner 
regions, and the other in which gas first flows into 
the central regions of a halo and is then trans- 
formed, once in place, into stars. In the first 
case, the stellar clumps loose their orbital energy 
through dynamical friction by heating the dark 
matter halo an d hence smoothing it s inner den- 
sity c usp fe.g.. | E1-Zant et aJjl20pil iBertin et al 



20031: iMa fc Bovlan-Kolchinll2004l ). In the second 
case, the gas radiates away its orbital and ther- 
mal energy and is responsible for the adiabatic 
contraction of the dark matter halo, i.e., for a 
central density profile ste eper with respect to the 
e.g. , iBlumenthal et al. 



initial distribution 



19861: 



Jesseit et ai]|2002t iGnedin et ail 12004) . In realis 
tic models, both processes are expected to occur 
and their relative contribution determines the ac- 
tual distribution of dark matter in the centre of 
a halo. For massive ellipticals, smo othed parti- 



cle hydrodynamical simulat ions (e.g., iNaab et al 
2007t Ijohansson et al" 2009) seem to favor a star 
formation picture which is characterized by an 
initial dissipative collapse of gas with quick "in 
-situ" star formation (at redshift > 2), followed 
by an extended phase of significant accretion of 
stellar material (at redshift < 3). Following the 
simulated time evolution of the stellar mass as- 
sembly in a massive galaxy, it is possible to notice 
that a remarkable fraction of the stars observed at 
the effective radius at redshift z ~ was not born 
there, but has been captured and accumulated 
in this region. The accretion of these stars, also 
called dry (gas-poor) minor mergers, is thought 
to be very important for the structural evolution 
of massive early-type galaxies and, in particular, 
at the origin of the late size, mass, and density 



Naabetal 



evolu tion of early-type galaxies (e.g. 
20091 ). 

The correlations observed in the first and third 
panels of Fig. |4]between the projected dark matter 
density and stellar mass and effective radius may 
be connected to the discussed mechanism of dissi- 
pationless accretion of groups of stars formed out- 
side the galaxy central regions. According to this 
picture, galaxies with a large number of captured 
satellites of stars should on average show light 
distribution that are more diffused, hence larger 
values of the effective radii, compared to galax- 
ies with less stellar satellites. Moreover, galax- 



ies with high luminous mass may be those that 
have been able to accrete more external stellar sys- 
tems. A large number of stellar clumps located in 
the inner regions of the galaxies with the largest 
values of effective radius and luminous mass can 
have smoothed the central cusp of the dark matter 
density profile appreciably. As a consequence, in 
the dissipationless scenario of galaxy evolution an 
anti-correlation between the values of the central 
(surface) dark matter density and effective radius 
and stellar mass is likely to be expected. 

Other processes which can contribute to deter- 
mine the amount and concentration of dark mat- 
ter in the central regions of galaxies are feed- 
back from supernovae and active galactic nuclei 
(AGN). The latter is thought to be particularly 
important in galaxies similar to those in our sam- 
ple, i. e., characterized by high velocity dispersion s 
(e-g-, 



Silk fc R.eeslll998t ICiotti fc Ostrikerll200ll) 



In detail, feedback from AGN is believed to be 
able to heat the cold gas in the central regions 
of ma ssive halos and quenc h star formation there 



Sprineel et all 120051: ISazonov et ail 120051: 



Ciotti fc Ostrikerll2007l )~ Despite the overall good 



agreement on this last point, it is still debated 
which mode of AGN activity is most important. 
A merger-induced "quasar mode" is thought to 
lead to an initial starburst followed by a suppres- 
sion of star formation and an expulsion of the gas 
from the galaxy center, once the super massive 
black holes become sufficiently massive. Given 
the extensively studied relation between the mass 
of a black hole and the ste llar velocity disper- 
sion o f its host galaxy 



2000t iGebhardt et al.l 120001 ) . the effect of this 
mechanism on the dark matter distribution should 
be more pronounced in galaxies which exhibit 
large values of stellar velocity dispersion. Further- 
more, the findings of several other studies (e.g. , 
Richstone et~allll998t iKormendv fc Bender! [2009) 



Ferrarese fc Merritt 



that have focused on the orbital decay of binary 
massive black holes arising from dissipationless 
mergers of gas-poor galaxies seem to predict re- 
sults that point to the same direction. In fact, this 
process may be at the origin of the formation of 
the observed cores in the surface brightness distri- 
bution of high-luminosity early-type galaxies. The 
scouring of a core in the stellar mass distribution 
of a galaxy through the ejection of stars from the 
galaxy center by the orbiting of the black hole bi- 
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nary is expected to depend on the mass of the ul- 
timately coalesced black hole, hence, on the value 
of the stellar velocity dispersion of the remnant 
galaxy. 

The second panels of Figs. [3] and [4] do in- 
deed show that galaxies with higher central stel- 
lar velocity dispersions present larger fractions and 
higher (surface) mass densities of dark matter pro- 
jected within the effective radius. We then specu- 
late that AGN feedback and / or massive black hole 
merging may contribute to the origin of these cor- 
relations. 

We notice that the measurements obtained in 
this work offer the interesting opportunity to test 
the predictions of numerical simulations and to 
quantify the relevance of the different physical pro- 
cesses that are believed to play a role in galaxy 
formation and evolution. 

6. Conclusions 

In this paper we have selected a sample of 
nearly 1.7 x 10 5 massive early- type galaxies at 
redshift smaller than 0.33 from the SDSS DR7 to 
study their projected dark over total mass fraction 
and surface dark matter density within the effec- 
tive radius. We have obtained from the MPA/JHU 
galaxy catalog the stellar mass values, that were 
estimated by fitting the SDSS multi-band photom- 
etry with a large grid of composite stellar popula- 
tion models and adopting a Chabrier stellar IMF. 
By generalizing the results derived from gravita- 
tional lensing and stellar dynamics analyses on a 
subsample of lens galaxies, we have modeled the 
total mass distribution of the galaxies in the sam- 
ple with a one-component isothermal model with 
effective velocity dispersion equal to the central 
stellar velocity dispersion and, hence, measured 
the galaxy total mass values inside their effec- 
tive radii. Then, we have combined the stellar 
and total mass measurements to determine the 
amount and concentration of dark matter within 
the galaxy half-light radii. Our main results can 
be summarized in the following points: 

— The projected total mass values are to a very 
good approximation linearly proportional to 
the stellar ones. This implies an almost con- 
stant fraction of dark over total mass con- 
tained within the galaxy effective radii. In 
detail, the dark matter component accounts 



on average for approximately 60% of the 
two-dimensional total masses. 

— The average logarithmic value of the surface 
dark matter density (in units of M@ kpc -2 ) 
inside R e is 9.1. The galaxies with the 
largest fractions of dark over total mass 
have the largest surface dark matter densi- 
ties within the effective radii. 

— If the stellar IMF and the dark matter den- 
sity profile do not show significant variations 
with the total stellar mass of the sample 
galaxies, the tilt of the FP of massive early- 
type galaxies cannot be primarily ascribed to 
variations in their central dark matter con- 
tent. 

— The observed correlations of the values of 
/d(< Re) and Sd,_r s with those of M*, cr , 
R e , and £*,.K e provide positive evidence on 
the importance of dry minor mergers, AGN 
feedback, and/or binary black hole coales- 
cence in driving the structural evolution of 
massive early-type galaxies. 

We remark that a plausible and viable solution 
to a better understanding of the results presented 
in this study awaits for cosmological simulations 
that are able to model the physics of baryons at 
a more accurate level than reached so far. The 
reproduction through simulation of these observa- 
tional findings would provide an invaluable step 
forward in our insight into the main mechanisms 
that rule the initial formation and the subsequent 
evolution of massive early- type galaxies. 
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